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Si 

a,: 
cr 



Mankei Tsang£ 

Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 02139, USA 
(Dated: March 6, 2008) 

bb" 

The ultimate electric and magnetic energy densities that can be attained by ban- 

o" 

£N) ■ dlimited electromagnetic pulses in free space are calculated using an ab initio quan- 

'_ 

tized treatment, and the quantum states of electromagnetic fields that achieve the 



ultimate energy densities are derived. The ultimate energy densities also provide 
an experimentally accessible metric for the degree of localization of polychromatic 
photons. 
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Ultrafast lasers have become an indispensible tool in a wide spectrum of science, including 
2 ] , metrology B , laser fusion [4] , biological imaging B , biological surgery 



>• 1 nonlinear optics [l|. 



7|. A key to the success of such lasers is the extremely high peak energy 



t*"- . |6|, and chemistry 

r-- ■ 

density that they can achieve, as the moderate energy of each laser pulse can be confined 

■ 

within femtoseconds or even attoseconds and focused to a micron-sized area. The high 

OO ! 

energy density strongly enhances light-matter interactions, and is especially crucial to the 



study of relativistic nonlinear optics [2|. Given the importance of an ultrahigh optical energy 
density in a broad range of applications, the limit to which one can focus a broadband optical 
pulse in three spatiotemporal dimensions and maximize the energy density is a fundamental 
problem. 

Localization of electromagnetic pulses has been investigated both in the classical regime 
js] and the quantum single-photon regime B; [lo|. While these seminal studies have con- 
tributed to our fundamental understanding of electromagnetic energy localization, all of 
their predictions have not yet been experimentally verified because of the difficulty in im- 



plementing their proposed electromagnetic pulse so 



contains arbitrarily high frequency components 



utions. Most require a spectrum that 



and can be exceedingly difficult to 



realize due to the finite laser gain bandwidth or finite transparency range of optical compo- 
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nents. The use of a spontaneously emitting atom proposed in Ref. [lOj couples the properties 
of the emitted photon to those of the atom and does not seem to be generalizable to other 
situations. Moreover, the quantum studies 0, Q are mainly concerned with the decay of 
the energy density far away from the center of localization at an instant of time for one 
photon, and thus are not immediately relevant to the more practical problem of maximizing 
the energy density at the center of localization for a large number of bandlimited photons. 

In this Letter, the ultimate electric and magnetic energy densities that can be attained by 
bandlimited electromagnetic pulses in free space are calculated using an ab initio quantized 
treatment, and the quantum states that achieve the ultimate densities are derived. By 
taking into account all degrees of freedom of electromagnetic fields and explicitly limiting 
the bandwidth of the pulses, our result overcomes all the shortcomings of the aforementioned 
studies and is more applicable to experimental situations. Measuring the energy densities 
at a fixed point is also considerably easier experimentally than measuring the decay of the 
energy density at an instant of time, so the maximum achievable energy densities can be used 
as an alternative and more accessible metric for the degree of localization of polychromatic 
photons. Most importantly, the ultimate energy densities impose a fundamental limit to 
which a bandlimited optical pulse can be focused spatially and temporally, so the presented 
result should prove useful for designing ultrafast optics experiments. 

The procedure of calculating the maximum energy densities is similar to the one used to 



calculate the multiphoton absorption rate limit for monochromatic light in Ref. , except 
that here we generalize the procedure to polychromatic light, such that all degrees of freedom 
are taken into account and the treatment can be considered ab initio. We also calculate the 
maximum magnetic energy density and the corresponding quantum state, as the magnetic 
field can also play a significant role in relativistic nonlinear optics Q]. 

We first derive the ultimate electric energy density, since it is more important for most 
applications. Consider the quantized electric field operator in free space 121 ]: 

E(r, t) = j^Jz E / d * k (^) V2 £ ( fc > + H.c, (1) 

where a denotes the two transverse polarizations, e(k, a) is the unit electric-field polarization 
vector, uj = ck = c(k 2 + k 2 + k 2 ) 1 / 2 is the frequency, a(k,a) is the annihilation operator 
satisfying the commutation relation [a(k, a), a^(k', a')] = 5 3 (k — k')5 aa >, and H.c. is the 
Hermitian conjugate. To impose a limit on the bandwidth, it is necessary to describe 
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the optical modes in terms of the frequency variable. This can be done by changing the 
momentum-space coordinates (k x ,k y ,k z ) to normalized spherical coordinates (fl, 0,0): 

oj = u Q Q, k x = koQ sin 9 cos 0, 

k y = k Q sin 9 sin 0, k z = koQ cos 9, 

dk x dk y dk z = dttd9d(j) {k^ 2 sin 9) , a(k, a) = a(fi, 9, 0, a) (/^fi 2 sin 0) " 1/2 . (2) 

where a;o is a normalization frequency, fco = ujq/c = 2tc/\q, and the annihilation operator has 
been renormalized so that its commutator is [a(Q, 9, 0, a), a'(Q', 9', 0', a')] = S(Q — Q')5(6 — 
9')6(<f) — (p')5 aa i. The positive-frequency electric field becomes 



E i+ \r,t) =i 



(hi \ ^ 2 r°° r n 
M / dVt d9 d<j)(Q 3 sm9) 1/2 e{9,<p,cx)d{Q,9,(f),a)e i 
2eoA / Jo Jo Jo 

(3) 
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In the continuous Fock space representation 
given by 

Y 



12j, the iV-photon momentum eigenstate is 



|^i,0i,0i,0i, • • ■ ,^N,9 N ,(p N ,a N ) = --L= Y\_d ] (Q n , 9 n , (p n , a n )\0) , (4) 
and the identity operator is 

oo 

i = J2\N)(N\, (5) 

N=0 

\N)(N\ = J dMideidfa . . . dtt N d9 N d0 N 

Cl,...,<TjV 

X |^1, 01, 01, 01, • • • , 0AT, 07V, 0jv)(^l, 01, 01, 01, • • • , ^N, 0JV, 0/V, &n\- (6) 

An arbitrary quantum state of electromagnetic fields can thus be expressed as 

oo 

\^) = Y^C N \N), C N =(N\V), (7) 

N=0 

and a Fock state as 

\N) — J dVLid9id<pi ... dQ N d9 N d(fi N & N (Qi, 9i, 0i, cr 1 , Q N , 9 N , (fi N , a N ) 

(T1,...,(TJV 

X |fil, 01, 01, 01, ... , fi/V, 07V, 07V, 07v), (8) 



where 



$jv(^l, 9l, 01, (Ti, . . . , Q, N , 6 N , <j) N , a N ) = (fii, 6>i, 01, (Ji, . . . , fijv, 0/v, 07V, CTjv|./V) (9) 

is the iV-photon momentum-space probability amplitude, which must satisfy the normaliza- 
tion condition: 

J dVtidOidcj)! . . . dVL N d9 N d(fi N \^ N {VL u 6 u ^ u a u . . . , Sl^v, 6W, 0jv, 0"jv)| 2 = !, (10) 

<ri,...,a N 

and the symmetrization condition: 

• • , ^n, @ni 0ri, &ni • • • , ^m, 9 m , 0m; @m, • • • ) 

= $at(. . . , Q m , 9 m , 0m, <T m , ■ ■ ■ , ^n, 9 n , 0„, <t„, . . . ) for any n and m. (11) 

To impose a limited bandwidth (a < Q < (3) on the electromagnetic fields, we require the 
probability of photons existing outside the bandwidth to vanish: 

$at(Qi, 0i, 0i, <7i, . . . , ^at, 6»at, 0_/v, ctat) = for any Q„ < a or Q n > (3. (12) 

With the theoretical framework put forth, we now proceed to calculate a bound on the 
electric energy density (minus the zero-point energy density) given by 



EE 



e E^ ■ E M 



(13) 



A bound on the electric energy density is equivalent to a bound on the energy density for 
one component of the electric field: 



e \ 2 



( p .|)':) = ( ,((p.* H )( p .J W )), 



(14) 



where p is an arbitrary real unit vector, because U e and U' e are equivalent if we choose p to 
be parallel to the electric field. In terms of the momentum-space representation, 

U 'e = ^ J2 \° N \ 2N J dtt 2 d9 2 d(j) 2 . . . dn N d9 N d<j) N 



x 



^=0 a 2 ,—,(JN 

n/3 P7V r2w 



pp pit pZir 

dQ d9 d<f> i (ft 3 sin 9) 1/2 p -e{9,(f),a)e 

„ Ja JO JO L 
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X ®n(Q, 9, 0, (J, Q 2 , 02, 02, 02, • • • , ^N, 9 N , 07V, <7n) 



(15) 
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l N Clu n CMp n 



where the symmetric property of <&jv given by Eq. (fTTj) is used. By virtue of the Schwarz's 
inequality and the normalization condition given by Eq. (fT(J|) . 

U' e < ^ jr \C N \ 2 N I dn x dB x d(j)x . . . d{l N d9? 

n=o cr l ,...,a N 

x |$jv(^i, 0i,0i, en, ■ ■ . ,Q N ,9 N ,4> N ,a N )\' 

r(3 pit / t 2n 

xJ^ dQ de d(j) (tt 3 sin 9) \p ■ i 
_ J a Jo Jo 
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where (N) = J2 N |Cjv | 2 ^ is the average photon number. As expected, the bound on U' e 
does not depend on p, and is therefore also applicable to the total electric energy density. 
Defining the actual lower and upper frequencies as uj\ = au and u 2 = /3u> , respectively, 
and the corresponding wavelengths as Ai,2 = 2'kc/uji j 2, we obtain the central result of this 
Letter: 

This simple expression agrees with the intuition that the ultimate energy density is limited 
by the maximum energy of photons ((N) U0J2) divided by the smallest volume that the 
photons can occupy (Af). 

In the limit of a small bandwidth compared to the center frequency, we can let Au = 
0J2 — lui, u = (u\ + u)z)/2, /S.U <^ Uo, and obtain 

!fE.E:)< 2 - { N } ^ (18) 

which is a bound on the peak intensity in the slowly-varying envelope regime, and again 
agrees with the intuition that the highest intensity is achieved when the mean energy of 
the photons is focused to their minimum pulse width (27r/Aa;) and beam size (Ajj). This 
approximate bound also agrees with that derived in Ref. ll|, where the monochromatic 
approximation is made at the beginning. Beyond the slowly-varying envelope regime, the 
exact bound given by Eq. (1171) depends on the upper frequency to the fourth power, un- 
derlying the importance of high-frequency components in maximizing the energy density, as 
they have a higher energy as well as a smaller localization volume. 

The use of the Schwarz's inequality is reminiscent of the matched filter concept in com- 
munication theory [13]. In Eq. ffl5l) . the A^-photon amplitude can be regarded as the input 
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signal, and the expression in square brackets can be regarded as a filter transfer function in 
the measurement of the electric field. An iV-photon amplitude that achieves the Schwarz 
upper bound is one that is linearly dependent on the square-bracketed expression, or in 
other words, when the input signal matches the filter. Assuming a factorizable <3?at, the 
following iV-photon amplitude that achieves the ultimate electric energy density can then 
be obtained: 

N 



71=1 

-iC- 1 ' 2 ^ sinOy^p ■ e*(6, 0, a ) e - ifc n>+^o for a < ^ n < p, 
otherwise, 



C = y (/? 4 - a 4 ) . (19) 

This state produces the maximum electric energy density at (r ,t ), with the electric field 
at (r ,t ) polarized along p. 

To apply the above result to the classical regime, let 

Cn = e -W/^. (20) 



Together with a factorizable $at in Eq. ( TL9l) . the quantum state becomes a coherent state in 



the continuous mode representation 



12j, and the mean electric field is then given by 



Jo 



x (fi 3 sin^ 1/2 £(^,0, ( r)/ e (fi,^0,cT)e ifc ™' + H.c., (21) 

which, incidentally, must be an exact solution of the Maxwell equations. The Fourier trans- 
form of E(r,t) is proportional to uj 3 , and the classical power spectrum is then proportional 
to cu 6 within the allowed frequency band. 
Consider now the magnetic field operator: 

B(r, t) = J d 3 k V2 n x s(k, a)a(k, a)e^~^ + H.c, (22) 

where n = k/k. While the total magnetic energy must be the same as the total electric 
energy for photons in free space, it is not difficult to show that the magnetic energy density at 
(tq, to) is zero where the electric energy density is maximum for the state given by Eqs. f|T9l) . 
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To maximize the magnetic energy density instead, we can simply apply the same procedure 
as above to the magnetic energy density, which turns out to obey the same bound as the 
electric one: 

k^^l^wU-Af)' (23) 

The quantum state with the ultimate magnetic energy density is also similar to the electric 

case, 

N 
n=l 

fb = \ (24) 
I otherwise, 

where m is the unit vector of the magnetic field at (?"o,£o)- 

The ultimate electric and magnetic energy densities may be challenging to achieve ex- 
perimentally, as they require a power spectrum proportional to u 6 within the allowed band, 
spatial focusing in all directions with a specific angular spectrum, and polarization control. 
That said, the results set forth impose fundamental limits to which the energy densities 
can reach regardless of technological advances in the control of electromagnetic fields, and 
therefore should prove useful for designing ultrafast optics experiments. 
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